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Abstract. We generalize a new class of cluster type mutations for which ex- 
change transformations are given by reciprocal polynomials. In the case of 
second-order polynomials of the form x + 2cos7r/n + x _1 these transforma- 
tions are related to triangulations of Riemann surfaces of arbitrary genus with 
at least one hole/puncture and with an arbitrary number of orbifold points of 
arbitrary integer orders n . We propose the dual graph description of the corre- 
sponding Teichmiillcr spaces, construct the Poisson algebra of the Tcichmiillcr 
space coordinates, propose the combinatorial description of the corresponding 
geodesic functions and find the mapping class group transformations. 

1. Introduction 

Since their appearance, cluster variables |14| find applications in geometry. An 
important example of the cluster variables in the case where exchange polynomial 
is of the second-order is provided by A-lengths [23], [22] of curves that make ideal 
triangle partitions of fundamental domains of Riemann surfaces with punctures. 
These systems were generalized in [5], [S] to the case of Riemann surfaces with 
holes. At the same time, a combinatorial description of geodesic functions in terms 
of the dual variables, the shear coordinates, as well as their quantization, was 
developed in [5]. Amazingly enough, transition from punctures to holes does not 
effectively change the corresponding cluster algebra. Generalizations of Tcichmiillcr 
spaces of Riemann surfaces to the case of bordered Riemann surfaces [20] or ciliated 
Riemann surfaces [lOj were constructed. The corresponding cluster algebras were 
developed in [19j [11] [12] , whereas the geometrical pattern underlying the bordered 
Riemann surfaces was identified with that of Riemann surfaces with Z^-orbifold 
points in [2], [3], where the corresponding mutations (flips) in terms of the shear 
coordinates were constructed. These flips preserve the sets of geodesic functions; 
on the level of cluster variables the corresponding transformations were considered 
in [13] : the mutations of cluster variables were there again given by the standard 
two-term relations. In [6], the description of Tcichmiillcr spaces of Riemann surfaces 
with holes and with orbifold points of order two and three was given. 

It was, however, an ambiguity in identifying a boundary component with an 
orbifold point; in the present paper, we demonstrate that this ambiguity is due to 
choice of the type of an orbifold point corresponding to the given boundary com- 
ponent; we therefore provide the combinatorial description of Riemann surfaces 
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with holes and with orbifold points of arbitrary orders. We show that mutations 
for orbifold points of order greater than two are given by three-, not two-term 
transformations determined by a second-order reciprocal polynomial. We prove the 
Laurent phenomenon and positivity property for these transformations. The posi- 
tive coefficients of Laurent polynomials however become not necessarily integers in 
the presence of orbifold points of order greater than three. On the shear-coordinate 
side, we define the complete set of real- valued coordinates, construct all the geo- 
desic functions for such surfaces, all the mapping-class-group transformations, and 
prove the regularity condition, that is, that all elements of the corresponding Fuch- 
sian group are hyperbolic or parabolic ones except precisely elements conjugate to 
passing around orbifold points. We therefore have a regular (up to exactly the 
indicated orbifold points) Riemann surface with holes for any choice of the intro- 
duced real coordinates, and vice versa; these coordinates parameterize therefore the 
corresponding Tcichmiillcr spaces of Riemann surfaces with holes and with orbifold 
points. 

The mutations of orbifold triangulations are particular cases of more general 
construction. In [18) some generalization of cluster transformations were described 
that preserve Poisson bracket and have some universal properties. Until recently 
no applications of these transformations were found. We were happy to see man- 
ifestation of our formulas in mutations of orbifold triangulations. Motivated by 
that we propose a new algebraic construction of generalized cluster algebras with 
mutations given by reciprocal polynomials. Using the tools of the standard cluster 
algebra [15j , |16j , we prove that the Laurent phenomenon holds true in this case as 
well. For algebras of order greater than two, we do not know whether the positivity 
property holds in general; it however holds in all tested examples, so we formulate 
it as a conjecture. We also prove that generalized cluster algebras of finite type 
satisfy the same Cartan-Killing classification as the standard cluster algebras. 
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2. Generalized cluster algebra 

We briefly remind the definition of cluster algebra. 

An integer nx n matrix B is called skew-symmetrizable if there exists an integer 
diagonal n x n matrix D = diag(d\, . . . , d n ), such that the product BD is a skew- 
symmetric matrix, i.e., bijdj — —bjidi. 

Let P be a semifield equipped with commutative multiplication • and addition 
0. We assume that the multiplicative group of P is a free abelian group. P is 
a coefficient group of cluster algebra. ZP is the integer group ring, J 7 is a field 
of rational functions in n independent variables with coefficients in the field of 
fractions of ZP. T is called an ambient field. 
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Definition 2.1. A seed is a triple (x, p,B), where 

• P = (Px)xex, a 2n-tuple of elements of P is a coefficient tuple of cluster x; 

• x = {sci, . . . , x n } is a collection of algebraically independent rational func- 
tions of n variables which generates J- over the field of fractions of ZP; 

• B is a skew-symmetrizable exchange matrix. 

The part x of seed (x, p, B) is called cluster, elements Xi G x are called cluster 
variables, and B is called exchange matrix. 

Definition 2.2 (seed mutation). For any k, 1 < k < n we define t/te mutation of 
seed (x, p, f?) in direction k as a new seed (x', p', B 1 ) in the following way: 



6ij , if i = fc or j 

b + \b ik \b kj +b ik \b kj \ otherwise . 



( 2 - 2 ) <=\ PtU bik >0 ->+P k U 



if z A/ j 

X ■ J 

— — , otherwise. 



(2-3) V't - Vl 

(2.4) for^fc p'+jp'r = { f^H^' ? J H ^ 

We write (x',p',B') = ^ fc ((x, p, B)). Notice that ^ k (fj, k ((x, p, S))) = (x, p, £?). 
We say that two seeds are mutation- equivalent if one is obtained from the other by 
a sequence of seed mutations. Similarly we say that two clusters or two exchange 
matrices are mutation- equivalent. 

For any skew-symmetrizable matrix B we define initial seed 

(x,p,S) = ({xi,...,x n l{pf,...,p±lB), 

where B is the initial exchange matrix, x = {x\, . . . , x n } is the initial cluster, 
p = {pf, ■ ■ ■ ,Pn} is th e initial coefficient tuple. 

Cluster algebra A(B) associated with the skew-symmetrizable n x n matrix B 
is a subalgebra of Q(x\, . . . , x n ) generated by all cluster variables of the clusters 
mutation-equivalent to the initial seed (x,B). 

Cluster algebra A(B) is called of finite type if it contains only finitely many 
cluster variables. In other words, all clusters mutation-equivalent to initial cluster 
contain totally only finitely many distinct cluster variables. 

Two most important properties of cluster algebra are Laurent phenomenon [15] 
and finite type classification [TB] . More exactly, Lurent phenomenon states that any 
cluster variable is expressed as a Laurent polynomial in terms of the initial cluster. 
The remarkable finite type classification claims that cluster algebras of finite type 
are in one-to-one correspondence with the Dynkin diagrams of finite type. 

2.1. Generalized cluster transformations. Now we introduce more general 
cluster transformations. 

Assume that B is a skew-symmetrizable integer matrix such that all elements in 
its kth row are divisible by d k . Define f3 k j = b k j/d k - 

Lemma 2.3. Let B' = fJ-i(B) is obtained from B by mutation in direction I. Then, 
all entries b' k j of kth row of B' are divisible by d k . 
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Proof. The statement follows immediately from matrix mutation 12. II □ 

Let us fix dk for all k from 1 to n and we assume that all elements bkj of kth 
row of integer skew-symmetrizable matrix B are divisible by dk ■ 

For a collection pi = (pj ; o, • • • , £»;d 4 ) we define the polynomial 9i[pi)(u, v) = 
~52iLoPi;£ ulvdi ~ l ^ e a polynomial of degree di. The corresponding inhomogeneous 
polynomial we denote by Pi[pi] = Pi[Pi](t) — Y^eLoPi;^- Note that, 6i[pi](u, v) = 
u di Pi[Pi}{v/u). 

We define a generalized seed q of a generalized cluster algebra as a triple 
q = (x(g), p(q), B(q)), where x = (.Ti(g), . . . , x n {q)) is a n-tuple of cluster variables 
in seed q, p(q) = (pi (<?), . . . , p n (q)) is n-tuple of coefficient collections Pi(q), Pi(q) = 
(pi-o{q), ■ ■ ■ ,Pi-,di{q)) is a di + 1-tuplc of coefficients of 6i[q], and, finally, B(q) is an 
exchange n x n matrix. 

Generalized cluster mutations are described by the following formulas: 
Exchange matrix is mutated in direction k by the rule 12.11 Introduce Uj ; >o = 
Y[p e>0 Xg e , Uj ;< o = e <o ' Mutation of cluster variables is given by the rule 
({x$,{pi} 1 B>) = p k ({x i i{p l },B): 

, , ,_j Xi, ifi^k; 

{ °> Xi ~ 1 fl *("^>°-"*;<o) otherwise. 

Coefficients mutate by the following generalized rule: 
(2-6) p' k . e = Pk-d k -£ 

(2.7) for*^ P >./ P >, r = { { f h ^;^ /P r ?MS 5 

J ' t {Pk;d k ) k% Pi;r/Pi;j, if Ofei < 0; 

Remark 2.4. Note that if we assume that coefficients of all 9i do not change 
under mutation then the corresponding inhomogeneous polynomial pi is reciprocal 
of degree d^ i.e., t di p{l/t) = p(t). 

Theorem 2.5. (Laurent property) Any generalized cluster variable is a Laurent 
polynomial in initial cluster variables X{. 

Proof. The proof uses the " caterpillar lemma" [15] . 

Lemma 2.6. Assume that a generalized exchange pattern on T n _ rn satisfies the 
following conditions: 

(1) For any edge the polynomial P does not depend on x k and is not divisible 
by any Xi, i G [n]. 

(2) Each exchange polynomial has nonnegative coefficients 

(3) For any three edges labeled by i, j, i 

i J i 

* P * Q * R * 

we have L ■ Q\ ■ P = R\ qq , where b is a negative integer, Qq — Q\ Xi =o, 

x i *~ 

and L is a Laurent monomial whose coefficient lies in A and is coprime 
with P. 

Then each element Xi(t) fori G [n], t G T n . m is a Laurent polynomial mxi(to), ■ • ■ ,x n (to) 
with coefficients in A. 
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By definition of generalized cluster transformation P = 8i[t], Q = 6j[t'}, R = 
9i[t"]. 

Note that parts (1) and (2) are evidently satisfied by generalized cluster muta- 
tions. It remains to proof part (3). note that if x% is not included into any monomial 
of Q then generalized mutation with labels i and j commute and the latter mutation 
is inverse to the former, namely, Xi(t"') = Xi(t). 

Consider the case when Xi enters a monomial of Q. For simplicity denote bji(t') 
by a. 

By our assumptions a ^ 0. Moreover, without loss of generality we can assume 
that a > 0, otherwise we replace B by —B. Since Q is determined by the homoge- 
neous polynomial of two variables where only one variable contains a positive power 
of Xi then Xi enters all monomial of Q but one. Hence, Qo = Q\ Xi =o is a monomial. 
Moreover, Q Q = p rfi II&. fc (t')>o x k ^* ^ Notc tnat P = ®i( u >o(t'), u <0 (t')). By the 

ba{t'), if 6,7(0 < 0; 

bu(t') + abji(t'), otherwise. 



mutation rule 12. 2\ 

b U (t" 



By the definition of generalized cluster transformation R = 0i(u>o(O> u<o(t"))- 

">o(g) 

«<o(?) ' 



For q = t' or t" introduce r q = Y[ t x^qf^ = H>°M 



Finally, 




-a/di 1 

Pj',0 T t 



It is enough to notice that 12.51 and 12.61 imply that pt;i(rt) = pt";i \ Tt"\ x , Qo 
Therefore, R\ q = P ■ L, where L is a Laurent monomial. □ 

Theorem 2.7. Generalized cluster algebras of finite type satisfy the same Cartan- 
Killing classification as the standard cluster algebras. 

Proof. The proof repeats the one of [TB]. The only differences make the proofs of 
the fact that the only finite type generalized cluster algebras of rank two correspond 
to A2, B2, G2 types. It is checked by direct computation similar to one in [TB]. Note 
first that in j42-case formulas for generalized cluster transformation coincide with 
formulas for the standard cluster transformation. In i?2-case the polynomial degrees 
of theta-polynomials arc two and one. Set the theta polynomials in the initial 
cluster 9i(u, v) = au 2 + buv + cv 2 , 62 (u, v) = pu + qv. Then, we immediately obtain 
(x,y) ^ ^ (x 1 ,y 1 ) ^ (x 2 ,yi) ^ (^2,2/2) ^ {x, y 2 ) ^ 

(x, y), where 

xi = [a + by + cy 2 )/x 7 

yi = (px + qa + bqy + cqy 2 )/xy, 

X2 = (a 2 q 2 + 2apqx + acq 2 y 2 + abq 2 y + bpqxy + p 2 x 2 )/xy 2 , 
2/2 = (qa + px)/y. 



(i 
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Similar computations lead to the 8-cycle in GVcase. Note that degrees of polyno- 
mials deg(0i) = 3, deg(#2) = 1- We set 6\{u, v) = au 3 + bu 2 v+cuv 2 +dv 3 , 02(u 7 v) = 
pu + qv. (x,y) ^ {m)(x lt y) ^ (xi,yi) ^> (^2,2/1) ^> (2=2,2/2) ^ 

(23,2/2) (x3,ys) ^> (x,ys) (x,y), where 
xi = (a + by + cy 2 + dy 3 )/x, 
2/1 = (px + aq + bqy + cqy 2 + dqy 3 )/xy, 

x 2 = (a 3 q 3 + 2a 2 cq 3 y 2 + 2a 2 dq 3 y 3 + 3a 2 pq 2 x + 2a 2 bq 3 y + 2abdq 3 y 4 + 3ap 2 qx 2 + 
Aabpq 2 xy + ac 2 q 3 y A + 2acdq 3 y 5 + ad 2 q 3 y 6 + 3acpq 2 xy 2 + 3adpq 2 xy 3 + ab 2 q 3 y 2 + 
2abcq 3 y 3 + bcpq 2 xy 3 + bdpq 2 xy 4 +p 3 x 3 + b 2 pq 2 xy 2 + 2bp 2 qx 2 y + p 2 cqx 2 y 2 )/x 2 y 3 , 

2/2 = {q 2 a 2 + abq 2 y + acq 2 y 2 + adq y + 2apqx + bpqxy + p 2 x 2 )/xy 2 , 

£3 = (a 2 dq 3 y 3 + a 2 cq 3 y 2 + acpq 2 xy 2 + a 2 bq 3 y + 2abpq 2 xy + bp 2 qx 2 y + a 3 q 3 + 
3a 2 pq 2 x + 3ap 2 qx 2 + p 3 x 3 )/xy 3 , 

2/3 = (px + aq)/y. 

□ 

We note that according to [T5] the generalized cluster transformation preserves 
presymplectic structure compatible with cluster algebra structure. Similarly, the 
secondary generalized cluster transformation preserves the compatible Poisson bracket. 

Theorem 2.8. (\18\ ) Poisson structure compatible with a cluster algebra is com- 
patible with the corresponding generalized cluster transformations. 

In the next section we describe generalized cluster structure associated with 
triangulated surfaces with orbifold points. 

3. TEICHMULLER SPACE OF SURFACES WITH HOLES AND ORBIFOLD POINTS OF 

ARBITRARY ORDER 

Wc now demonstrate how the above mutations with reciprocal polynomials of the 
second order appear in the description of Teichmullcr spaces of Riemann surfaces 
of arbitrary genus with nonzero number of holes (punctures) and with an arbitrary 
number r of orbifold points of arbitrary orders. 

3.1. The ideal triangle decompositions of orbifold Riemann surfaces and 
cluster variables. We now present the geometric pattern underlying the algebraic 
construction of cluster variables corresponding to orbifold Riemann surfaces. 

Particular cases of Riemann surfaces with orbifold points of order 2 and 3 are dis- 
cussed in 2, 4 , and [6] . For relation between skew-symmetrizable cluster algebras 
of finite mutation type and surfaces with orbifold points see |13j . 

We consider a regular Riemann surface S S)Sir of genus g with s > holes and 
with a number r > of orbifold points; orders of these points Pi, i = 1, . . . ,r are 
positive integers greater than one. 

We introduce the marking on the set of orbifold points splitting this set into 
nonintcrsecting subsets 5k (maybe empty), k = 1, . . . , s, J2k=i l^fcl = r - We then 
assign the subset Sk to the fcth hole for every k and introduce a cyclic ordering 
inside each subset 5k- 

To construct the generalization of the ideal triangle decomposition [51 [221, we nrs t 
remove from the surface all the hyperbolic domains of holes bounded by the cor- 
responding perimeter geodesic lines (with their closures, which are these perimeter 
lines). Second, we choose for each orbifold point from the set 5k a domain contain- 
ing this point and bounded by a geodesic curve whose both ends spiral to the fcth 
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FIGURE 1. The right part represents an example of a regular genus 
zero Riemann surface with a Z4-orbifold point marked by a star; the 
bounding geodesic line is the image of a side of an ideal equilateral 
square in the Poincare disc depicted in the left part. Both ends of the 
geodesic line spiral asymptotically to the closed geodesic that is the 
boundary of a hole. 



hole as shown in Fig. We remove from the remaining part of the surface all such 
domains. The remaining part of £ g s r admits splitting into ideal triangles; a copy 
of this splitting can be drawn as a connected ideal polygon in the Poincare disc; the 
sides of this polygon are of two sorts: those that are not pre-images of geodesies 
going around orbifold points must be pairwise identified; to a side that is a pre- 
image of the geodesies going around Z p orbifold point we attach (from outside) an 
equilateral ideal p-gone with the orbifold point situated at its geodesic center; this 
p-gone is the p-fold covering of the removed domain enclosing the orbifold point . 

An example of a fundamental domain of the Riemann surface £1,1,2 of genus one 
with one hole and with two Z3 orbifold points is in Fig. [5J 

We use the standard geometric correspondence between cluster variables and A- 
lengths: at each point at the absolute that is a vertex of an ideal triangle we set an 
horocycle; the A-length l a is then the geodesic length of the part of side a enclosed 
between two horocycles based at its endpoints. The correspondence reads: 

(3.8) a = e^'\ 

where £ a is understood as the signed distance between horocycles: £ a is negative 
when the corresponding horocycles overlap. Here a is the cluster variable associated 
with the edge. 

For an ideal quadrangle with the (cyclically enumerated) sides e^, i = 1, ... ,4 
and with diagonals d and d! we have the standard two-term cluster relation 

(3.9) dd! = ai<23 + a-ia^. 

We now consider mutations for cluster variables of the ideal p-gone corresponding 
to a Zp-orbifold point. We consider the pattern in the left-hand side of Fig. [3] and 
perform a sequence of mutations (|3.9|) to come to the pattern in the right-hand 



In the case of a Z2-orbifold point, this domain has zero area because the corresponding 
geodesies goes directly to the orbifold point, reflects at it, and repeats its path in the opposite 
direction. 
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FIGURE 2. The Poincare disc with depicted fundamental domain for 
the genus-one surface £1,1,2 with one hole and two Z3 orbifold points 
marked by *). Solid lines constitute the boundary of a fundamental 
domain and dashed lines are sides of the related ideal triangles. The 
sides with labels a and b are pairwise identified. Inscribed circles with 
centers at the Z3 orbifold points all have the radius | log 3. 




FIGURE 3. The mutation in the ideal p-gone (here p = 6), which is 
the p-fold covering of the domain around a Z p -orbifold point (marked 
by *). We perform standard 2-term mutations (|3.9p on the set of cluster 
variables to come from the pattern in the left-hand side to the one in 
the right-hand side. 



side. For a equilateral p-gone, the cluster variables for fc-diagonal are 

, 91 m sin(7rfc/p) 

(3.10) c k = c , (ci=c). 

sm(7r/p) 

The easy combinatorics then yields 

Lemma 3.1. X-lengths satisfy the following relation 

cc' = a 2 + 2 cos(-7r/ p)ab + b 2 . 

Remark 3.2. Relation in Lemma T3.ll is a generalized cluster mutation described 
in Sec. [5J Note that since we use only two-term transformations 13.91 to prove 
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Lemma I3.R the positivity property for generalized mutations of such form follows 
from the one for the standard mutations. 

To simplify the description, it is convenient to introduce the notion of petal 
surface. Petals are the domains containing orbifold points that were removed on 
the second step of constructing the ideal triangulation. To the fcth hole, we associate 
the bouquet of \5k\ petals (with no petals if Sk is empty), each petal carries, besides 
its cluster variable, the number uj p = 2cos(-7r/p). The mutation then occurs inside 
the corresponding ideal triangle (painted by a light color) , and we have three cases 
depending on whether the adjacent sides are petals themselves: 




I II 




In all the three cases above, the mutation law is given by Lemma \3. II In the first 
case, we transfer the cluster variable from one set Sk to another set (if the vertices 
V\ and V2 are distinct); these transformations may also change the cyclic ordering 
of orbifold points inside a set Sk- Note that the label u) p remains assigned to the 
transformed edge. 

3.2. Fat graph description for Riemann surfaces with holes and with Z p 
orbifold points. In this subsection and in the rest of the paper, we use the graphs 
dual to the above ideal triangle decompositions of Riemann surfaces. These graphs 
are especially useful when describing the Fuchsian groups A g s r of Riemann surfaces 
£ g , s ,r and the corresponding geodesic functions. 

Definition 3.3. Wc call a fat graph (a graph with the prescribed cyclic ordering 
of edges entering each vertex) T 9tS _ r a spine of the Riemann surface S S)Sir with 
g handles, s > holes, and r orbifold points of the corresponding orders pi, i = 
l,...,r, if 

(a) this graph can be embedded without self- intersections in £ 9 ,s. r ; 

(b) all vertices of T g , s ,r are three-valent except exactly r one-valent vertices 
(endpoints of "pending" edges) , which arc placed at the corresponding orb- 
ifold points; 

(c) for an orbifold point from the set Sk, the corresponding pending edge pro- 
trudes towards the interior of the face of the graph containing the fcth hole; 
the cyclic ordering of pending edges pointing towards the interior of this 
face coincide with that of orbifold points in the set Sk, 
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(d) upon cutting along all edges of r 9jS . r the Riemann surface £ g . s ,r splits into 
s polygons each containing exactly one hole and being simply connected 
upon gluing this hole. 

Edges of this graph are labeled by distinct integers a = 1, 2, . . . , 6g — 6 + 3s + 2r, 
and to each edge we set into the correspondence the real number Z a . 

The first homotopy groups 7Ti(£ ffiS)T .) and 7Ti(r giSjr ) coincide because each closed 
path in S g . Sj r can be homotopically transformed to a closed path in T g jS . r (taking 
into account paths that go around orbifold points) in a unique way. The standard 
statement in hyperbolic geometry is that conjugacy classes of elements of a Fuchsian 
group Ag^r are in the 1-1 correspondence with homotopy classes of closed paths 
in the Riemann surface £ s , s .r = H+/A fliSir . and that the "actual" length i 1 of a 
hyperbolic clement 7 G A fliSjr coincides with the minimum length of curves from 
the corresponding homotopy class; it is then the length of a unique closed geodesic 
line belonging to this class. 

When orbifold points arc present, the Fuchsian group contains besides hyper- 
bolic elements also elliptic elements corresponding to rotations about these orbifold 
points. The corresponding generators Fi, i = l,...,r, of the rotations through 
2tt/ pi are conjugates of the matrices 

(3.11) F i = U i F Pi Ur\ F p =( ^ _M, w = 2coS7r/p, p > 2. 



The real numbers Z a in Definition 13.31 are the h- lengths [22] : they are called 
the (Thurston) shear coordinates |23j.[Tj in the case of punctured Riemann surface 
(without boundary components). We preserve this notation and this term also 
in the case of orbifold surfaces. These coordinates are related to the cross-ratio 
relation for two adjacent ideal triangles constituting the ideal quadrangle with the 
respective vertices (in the cyclic order) a, b, c, d and diagonal bd. At the same 
time, they are related to the cluster variables corresponding to the sides of the 
corresponding quadrangle. We have 



(3.12) e 



(b - c)(d - a) _ aia 3 



(b — a)(d — c) 0204 



and we obtain the parameter Z a choosing {a, b, c, d} = { — 1, 0, e Za , 00}. 

For example, the fat graph corresponding to the pattern in Fig. [5] is depicted in 
Fig.|U 

3.3. The Fuchsian group A g s r and geodesic functions. We now describe 
combinatorially the conjugacy classes of the Fuchsian group A 9iS:r . Every time the 
path homeomorphic to a (closed) geodesic 7 goes through the edge with the label 
a we insert [5] the matrix of the Mobius transformation 



(3-13) X Za = 



-e z °/ 2 
e -z a /2 o 



We also have the "right" and "left" turn matrices to be set in proper places when 
a path makes the corresponding turns at three-valent vertices, 

(3.14, *_(_', I), L-tf-f. , \ 
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FIGURE 4. The fat graph corresponding to the ideal triangle partition 
of the fundamental domain in Fig. [3] The real numbers Z\ and Z2 are 
associated to the pending edges and A, B, and I2, Y3, and Y4 to the 
inner edges. 

New elements of the Fuchsian group correspond to rotations of geodesic lines 
when going around orbifold points indicated by star-vertices; for a Z p orbifold point 
we then insert the matrix F p (|3.11|) into the corresponding string of 2 x 2-matrices 
(when we go around the orbifold point counterclockwise as in Fig. Ufa)). When 
the order of the orbifold point is larger than two, we can go around it k times; we 
then have to insert the matrix (—l) k+1 Fp into the product of 2 x 2-matrices. For 
example, parts of geodesic functions in the three cases in Fig. [S]read 

(a) ...X x LX z F p X z LX Y 

(3.15) (b) ...XxLXzi-F^XzRXx..., 
(c) ...X Y RX Z (F^)X Z LX Y .... 

Note that F? = (— 1) P_1 E, so going around the Z p orbifold point p times merely 
corresponds to avoiding this orbifold point due to the simple equality (note that 
X 2 S = — E and L 2 = —R) 

XxLXzi^ — 1)^ ^ F^ X Z LX Y — XxLX^LXy — — XxL^Xy = XxRXy- 

(For the Z2 orbifold points this pattern was first proposed by Fock and Goncharov 
|10) : the graph morphisms were described in [3].) 

An element of a Fuchsian group has then the typical structure 

(3.16) P 7 = LX Yn RX Yn _ 1 ■ ■ ■ RX Y2 LX Zl {-l) k+1 F^X Zl RX Yll 

where Yi are variables of "internal" edges and Zj are those of pending edges. The 
corresponding geodesic function 

(3.17) G 7 = TrP 7 = 2cosh(£ 7 /2) 

is expressed via the actual length £ 1 of the closed geodesic on the Riemann surface. 
Remark 3.4. Note that the combinations 

p Y ( z- Yl2 ~e Y/2 \ 1 t y ( e- y li \ 

UJl V — I q £ Y/2 I ancl — I _ e -Y/2 e Y/2 J 
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FIGURE 5. Part of a graph with a pending edge. Its endpoint with the 
orbifold point is directed toward the interior of the boundary compo- 
nent this point is associated with. The variable Z corresponds to the 
respective pending edge. We present four typical examples of geodesies 
undergoing single (a), double (b), and triple (c) rotations at the Z p 
orbifold point. 



as well as products of any number of these matrices have the sign structure ( ^ 

so the trace of any of P 7 in the absence of orbifold points is a sum of exponentials 
with positive integer coefficients; this sum always include the terms e Yl l 2 ^ ^n/ 2 
and e~ Yl l 2 Y ™/ 2 being therefore always greater or equal two thus describing a 
hyperbolic or parabolic clement; the latter is possible only if Y% + • • ■ + Y n = 
and only if the turn matrices in (|3.16p are all either R or L, which corresponds to 
a path going along the boundary of a face; all such passes are homeomorphic to 
the hole boundaries, and the condition that the sum of Yi equals zero indicates the 
degeneration of a hole into a puncture. 

The group generated by elements (|3.11|) together with translations along A- and 
.B-cycles and around holes not necessarily produces a regular (metrizable) surface 
because its action is not necessarily discrete. We formulate the necessary and 
sufficient conditions for producing a regular surface in terms of graphs (see [3] for 
the Z2 orbifold point casc)0 

To formulate the regularity condition, we interpret passages around orbifold 
points as paths in the p-fold covering of the geodesic neighborhood in Fig. [TJ For 
this, we take the subgraph in Fig. [6] dual to the ideal triangle decomposition in 
Fig. [3] When splitting the equilateral p-gone into ideal triangles we break the p- 
fold symmetry, so now the shear coordinates Zi, i = 1, . . . , p, on the fat graph edges 
dual to the corresponding p-gone sides and Yj, j = 2, . . . ,p — 2, on the edges dual 
to the diagonals of the p-gone are different. 

We first identify the parameter Z in (|3.15[) to be 



z _ sin (t/p - 0/2) 



(3.18) 

sm(<p/2) 

We can then derive the explicit relations between the parameter Z in (|3.18[) and the 
variables Zi and Yj determined by the standard cross-ratio relations (|3.12j) . The 
vertices of the ideal p-gone are situated at the points e l27rfe / p , k = 0, . . . ,p — 1, and 



In what follows, wc call a Ricmann surface regular if it is locally a smooth constant-curvature 
surface everywhere except exactly r orbifold points. 
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FIGURE 6. The (tree-like) subgraph dual to the ideal triangle decom- 
position of the ideal equilateral p-gone in Fig. [3] All the variables Z a 
and Y/3 are determined by the cross-ratio relations in the correspond- 
ing ideal quadrangles based on the points from the sets {e 27rlk,lp } and 
{ e ^«=/p+#} j fe = 0,...,p-l. 



p copies of the vertex of an additional ideal triangle adjacent to the p-gone side are 
e it+i2wik/ P } k = 0,...,p-l. 

Using the cross-ratio relations to calculate Zi and Yj we find the exact relations 
between these variables and the variable Z given by (|3.18l) : 



e 



z l = c z sin(27r/p) 
sin(7r/p) ' 
c z sin(7r/p) 
sin(27r/p) ' 



v ? sinfffc — l)ir/p) 

(3.19) e z » = e z y ) ' P> , k = 2,...,p-l, 

sin(K7r/p) 

v sinf (k + l)ir/p) 

(3.20) e Y " = . 7 , k = 2,...,p-2. 

sm((fc — l)7r/pj 

The following 2 x 2-matrix equalities can be verified directly: 
(3.21) 

XzF p Xz = Xz±LXz 2 = Xz k LXy k LXz k+1 = Xz p _ 1 LXz p , k = 2, . . . ,p — 2. 
We then have the following lemma. 

Lemma 3.5. We have the following explicit 2 x 2-matrix relations for the shear 
variables in the equilateral p-gone in Fig. [S| given by h3.19)) and A3.20\) : 

XzFpXz = Xz 1 LXz 2 
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Xz(— Fp)Xz = Xz 1 RXy 2 LXz 3 

(3.22) : 

X z (-l) k - x F*X z = X Zl RX Y2 R-- RX Yk LXz k+1 , k = 2,...,p-2 
Xz(—1) P FP 1 Xz = Xz 1 RXy 2 R- • ■ RXy p _ 2 RXz p - 

The proof uses equalities from (|3 . 2 1 1) for constructing the longer chain using 
that XsXs = — E for any variable S and that L 2 = —R. For example, we obtain 
the r.h.s. of the second equality in (|3.22[) multiplying Xz 1 LXz 2 ■ Xz 2 LXy 2 LXz 3 
whereas the l.h.s. merely becomes XzF p Xz ■ XzF p Xz = Xz{—^)F p Xz- AH other 
equalities are obtained if we continue this chain of multiplications. 

Due to Lemma \'3. 51 all "rotations" about orbifold points Xz{— \) k+1 F p Xz are 
now presented as the standard products of matrices Xs (with real S) alternated 
with the matrices of left and right turns (|3.14p . which means that the conclusion of 
Remark l3.4l is valid in this case as well: as soon as in the original spine F g , s ,r all the 
parameters Z a of pending edges are real, all the geodesic functions constructed on a 
Ricmann surface with orbifold points are Laurent polynomials with positive integer 
coefficients of the "new" real variables z\ a ' p \ (where the superscripts a,p 

indicate that these variables are completely determined by the original variable Z a 
and the order p of the orbifold point) and "old" variables of "internal" edges of the 
spine r ffjS r . corresponding to usual partitions into ideal triangles. So, again, in the 
trace of every product of form (|3.16p we necessarily have the term 2 cosh(^^ =1 Xp), 
where the sum ranges all edges (new and internal ones) the corresponding path 
goes over, and we let X denote the variables on all these edges disregarding their 
origins. Every such trace is therefore a positive number greater or equal two, and 
the corresponding clement of the group will be either hyperbolic or parabolic (the 
latter is possible only for geodesies around holes and only if a hole reduces to a 
puncture). The only elliptic elements are precisely conjugates of F™. We therefore 
come to the theorem 

Theorem 3.6. We have a metrizable Riemann surface for any choice of real num- 
bers Z a associated to the edges of an original spine T g ^ Sjr . The converse statement 
is also true: for any metrizable Riemann surface £ 9 , s ,r we have a spine r gjSjr with 
real numbers associated to its edges such that the lengths of geodesies on £ giS , r are 
given by traces of products H3.16}) corresponding to paths in the spine. 

The proof of the second statement was performed in j4j for Z2-orbifold points. It 
is based on the (obvious) existence of the ideal triangle decomposition described in 
Sec. 13.11 for any metrizable Riemann surface and can be straightforwardly general- 
ized to the case of orbifold points of any type. We have therefore parameterized all 
possible regular surfaces in terms of the (6g — 6 + 3s + 2r)-tuple of real coordinates 
{Z a }. 

Corollary 3.7. The decorated Teichmuller space r of Riemann surfaces with 
holes and orbifold points is the space R 6 g-6+3s+2r Q j rea ^ p arame i ers on £/j e edges 
of a spine T g>s , r . 

4. Mapping class group transformations 

4.1. Poisson structure. One of the most attractive properties of the graph de- 
scription is a very simple Poisson algebra on the set of parameters Z a . The following 
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result is the straightforward generalization of the theorem formulated for surfaces 
without marked points in [8] and for surfaces with order-2 orbifold points in |10j 
(see also [5]). 

Theorem 4.1. In the coordinates Z a on any fixed spine corresponding to a surface 
with orbifold points, the Weil-Petersson bracket B WP is given by 



the cyclically (counterclockwise) ordered (ay = a\) edges incident to the vertex with 
the label a irrespectively on whether these edges are internal or pending edges of 
the graph. This bracket gives rise to the Goldman bracket on the space of geodesic 
length functions [17) . 

We identify the exchange matrix B with the matrix of the Poisson relations for 
the variables Z a . 

The center of this Poisson algebra is provided by the proposition. 

Proposition 4.2. The center of the Poisson algebra (|4.23[) is generated by ele- 
ments of the form ^2Z a , where the sum ranges all edges of Tg.s.r belonging to 
the same boundary component taken with multiplicities. This means, in particular, 
that each pending edge, irrespectively on the type of orbifold point it corresponds to, 
contributes twice to such sums. The dimension of this center is obviously s. 

for the proof in the general case see Appendix B of [7] . 

4.2. Flip morphisms of fat graphs. In this section, we present the complete list 
of mapping class group transformations that enable us to change numbers \5k\ of 
orbifold points associated with the fcth hole, change the cyclic ordering inside any of 
the sets 5k, flip any inner edge of the graph and, eventually, change the orientation 
of the geodesic spiraling to the hole perimeter (in the case where we have more than 
one hole)0 We can therefore establish a morphism between any two of the graphs 
belonging to the same class T g ^, r with the same (unordered) sets of orbifold point 
orders {pi] r i=1 - 

4.2.1. Whitehead moves on inner edges. Given a spine T of E and assuming that the 
edge a has distinct endpoints. we may produce another spine T a of E by contracting 
and expanding edge a of T, the edge labeled Z in Figure [71 This transformation is 
dual to the mutation (|3.9[) . We say that r a arises from T by a Whitehead move (or 
flip) along the edge a. A labeling of edges of the spine T implies a natural labeling 
of edges of the spine T a ; we then obtain a morphism between the spines T and T a . 

Proposition 4.3. [5J Setting 4>(Z) = log(l + e z ) and adopting the notation of 
Fig. for shear coordinates of nearby edges, the effect of a Whitehead move is 



(4.24) W z : {A,B,C,D,Z)^{A + 4>{Z),B-4>(-Z),C+<j>{Z),D-(j>{-Z),-Z) 



In the various cases where the edges are not distinct and identifying an edge with its 
shear coordinate in the obvious notation we have: if A = C, then A' = A + 2<fi(Z); 



if B = D, then B 1 = B - 2<f){-Z); if A = B (or C = D), then A' = A + Z (or 



where the sum ranges all the three-valent vertices of a graph and on are the labels of 




^These transformations are dual to mutations of cluster variables from Sec. [2] 
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FIGURE 7. Flip, or Whitehead move on the shear coordinates Z a . The 
outer edges can be pending, but the edge with respect to which the 
morphism is performed must be an internal edge. We also indicate the 
correspondences between geodesic paths under the flip. 



C' = C + Z);ifA = D (or B = C), then A' = A + Z (or B' = B + Z ). Any subset 
of edges A, B , C , and D can be pending edges of the graph. 

We have the lemma establishing the properties of invariance w.r.t. the flip mor- 
phisms [5]. 

Lemma 4.4. Transformation (|4.24[) preserves the traces of products over paths 
(|3.17p ( the geodesic functions ) and transformation (|4.24[) simultaneously preserves 
Poisson structure (|4.23[) on the shear coordinates. 

4.2.2. Whitehead moves on pending edges. Choosing other representatives of the 
orbifold points in the Poincare disc, we obtain different fundamental domains with 
different cyclic ordering of the (preimages) of the orbifold points Sj (i = 1, . . . , \5k\) 
possibly with transferring orbifold points from one set 5k to another set 5k' ■ 

Analogously to the mutation in Fig. [3l flipping the pending edge corresponds to 
choosing another fundamental domain, as shown in Fig. [5J We take there e z given 
by formula (|3.18[) and e Yl ' 2 and e Fl > 2 given by the standard cross-ratio relations, 
for example, 

^ (\ _ e i<f>^ ( £ ib _ e 27ri/pj ^ e -2iri / p+ib _ j\ ^g-27ri / p+i<f> _ gi<^ 
gifc^ g2iri/p^ ' ^g — 2-!ri/p+i<f> g—2iri/p-\-ib^ e'^) 

Lemma 4.5. The transformation in Fig. []| with e z given by i3.18\) has the form 

(4.25) {Y 1: Y 2 , Z} = {Y x - log(l + we' 2 + e- 2Z ),Y 2 + log(l + we z + e 2Z ), -Z} 

and is the morphism of the space 7~ g H sr . These morphisms preserve both Poisson 
structures (|4.23j) and the geodesic functions. In Fig. [3 any (or both) of Y -variables 
can be variables of pending edges (the transformation formula is insensitive to it). 

Proof. Verifying the preservation of Poisson relations (|4.23[) is simple, whereas 
for traces over paths we have four cases, and in each of these cases we have the 
following 2 x 2-matrix equalities to be verified directly: 

X Y2 LX z F<;{-l) k - 1 X z LX Yl = X Y2 RX 2 FP- k+1 (-iy- k X 2 RX Yi , 

X Yl RX z F^-l) k - 1 X z LX Yl = X Yi LX 2 FP- k+1 (-l)P- k X 2 RX Yi , 

X Y ,LX z Fl(-\f- x X z RX Y , = X Y RX- z Fl- h ^{-\y- k X- z LX % . 
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IT 




FIGURE 8. The transformation of dual variables (/i-lengths) 
{Yi, Y 2 , Z} {% Y 2 , Z} described by (|4/25)l with w = 2 cos(tt/p). 



Using flip morphisms in Fig. [5] and in formula (|4.24[) . we establish a morphism 
between any two algebras corresponding to surfaces of the same genus, same number 
of boundary components, and same numbers of Zp-orbifold points of each sort p; the 
distribution of latter into the boundary components as well as the cyclic ordering 
inside each of the boundary component can be arbitrary. 

It is a standard tool that if, after a series of morphisms, we come to a graph of 
the same combinatorial type as the initial one (disregarding labeling of edges but 
distinguishing between different orbifold types of pending vertices) , we associate a 
mapping class group operation to this morphism therefore passing from the groupoid 
of morphisms to the group of modular transformations. 

Remark 4.6. Another way of interpreting transformations f|4. 25[) is as follows. We 
can imitate the above flips/mutations by introducing a new "hole" with possibly 
imaginary perimeter P and considering the following chain of standard flips: 




Y- log(l-e- z )- log(l+e- p - z ) 



In this pattern, the variable Z is not completely defined by the ideal triangle 
decomposition: there is an arbitrariness in choosing the "mirror" triangle (we must 
just choose it consistently). So, in the above sequence of flips, it is useful to shift 
the variable and consider 

_ P 

z = z + T 
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The transformation for the variables X and Y then just becomes (|4.25[) . 

X + log[(l + eZ- p / 2 )(l + e^+ p ' 2 )} 



X 
Y 
~Z 



y-io g [(i 
-z 



,-Z+P/2 



)(1 



14 

r-iog[i 



log[l +uj p e z + e 2Z ] 



-2Z~\ 



where ui. 



-P/2 



Since P is not affected by the above sequence of flips, we can merely erase 
the corresponding loop and present it as flipping the pending edge, which carries 
besides the cluster variable Z also the "index" variable w p , which is not affected by 
mutations and which is equal to 2cos(7r/p) in the geometric case. 

4.2.3. Changing the spiraling direction. The last mapping class group transforma- 
tion changes the sign of the hole perimeter: 

P -P 



(4.26) 





Lemma 4.7. Transformation 
geodesic functions. 



preserves the Poisson brackets and the set of 



Proof. The preservation of the Poisson bracket is obvious because the variable 
P Poisson commutes with all other variables, whereas the preservation of geodesic 
functions follows from two matrix equalities: 

Xy LX pLXy — Xy -]-P LX _ pLXy_L_p, 
XyRXpRXy = Xy+pRX—pRXy+p. 

We can therefore enlarge the mapping class group of 7~ g H s r by adding symme- 
tries between sheets of the 2 s -ramified covering of the "genuine" (nondecorated) 
Tcichmullcr space T g . s , r - 

The geometrical meaning of this transformation is clear: we change the direction 
of spiraling to the hole perimeter line for all lines of the ideal triangle decomposition 
that spiral to a given hole like in Fig. Q] 

We can summarize as follows. 

Theorem 4.8. The whole mapping class group o/S ffiS r is generated by morphisms 
described by Lemmas\4-.4\ [7^1 and\4-.7\ 



5. Conjectures 

Lemma |3 . f I shows that generalized transformation a 2 + 2cos(ir/ p)ab + b 2 appears 
as a flip in the presence of an orbifold point of order p. The generalized cluster 
algebra constructed in this way is a subalgebra of a bigger standard cluster alge- 
bra (maybe of infinite rank) associated with triangulated surface while generalized 
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exchange relation are sequences of standard mutations. Note that the positivity of 
Laurent polynomials for cluster algebras associated with bordered surfaces is known 
by |21j . This implies the positivity of Laurent polynomials in generalized cluster 
algebra associated with triangulations of the surface with arbitrary orbifold points. 
We formulate the following conjecture. 

Conjecture 5.1. If p is a reciprocal polynomial with positive coefficients then any 
cluster variable of a generalized cluster algebra is expressed as a positive Laurent 
polynomial in the initial cluster. 

We checked by direct inspection that the statement holds for finite type rank 2 
cluster algebras. 

The example above leads to a natural question: 

Question: Is any generalized cluster algebra a subalgebra of some standard 
cluster algebra? 
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